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INTRODUCTION 


One  of  the  essential  properties  of  a  propellant  is  its  burning  rate  which, 
unfortunately,  has  been  poorly  defined.  Conceptually,  there  are  various  burning 
rates. 


The  strand  burner  gives  a  bulk  burning  rate  at  constant  pressure  but  it  may 
not  reflect  the  rate  at  which  the  surface  of  individual  propellant  grains  re¬ 
cedes,  nor  does  it  consider  any  interactive  effects  between  grains  or  dynamic 
effects  due  to  the  variation  of  pressure. 

A  burning  rate  determined  from  a  closed  chamber  firing  should  more  closely 
reproduce  the  conditions  actually  encountered  in  a  gun;  however,  current  methods 
of  evaluating  burning  rates  from  the  pressure  time  trace  do  not  yield  the  best 
parameters  for  use  with  ballistic  codes.  The  same  computer  ballistic  code  should 
be  used  to  evaluate  the  bomb  data  and  to  perform  the  ballistic  calculation;  that 
is,  a  closed  bomb  should  be  viewed  as  a  gun  with  an  infinite  shot  start  pres¬ 
sure.  In  this  way,  the  burning  rate  parameters  obtained  from  the  bomb  will  at 
least  partially  compensate  for  approximations  made  in  the  ballistic  code  and  thus 
compensate  for  the  same  approximations  made  when  calculating  the  muzzle  velocity. 

Finally  the  overall  burning  rate  is  determined  by  the  slowest  physical  or 
chemical  step.  The  faster  processes  cannot  proceed  until  this  slower  step  has 
been  completed.  The  specific  step  cannot  be  identified  at  this  time  but  it  is 
probably  in  the  gas  phase.  Otherwise,  there  would  not  be  evidence  of  an  exponen¬ 
tial  dependence  of  pressure  on  the  burning  rate. 

In  view  of  the  importance  of  the  burning  rate  in  propellants  and  the  fact 
that  there  are  several  phenomenological  definitions  of  this  quantity,  an  effort 
should  be  expended  to  clarify  the  conditions  and  to  explain  the  discrepancies 
found  in  the  different  methods  of  measuring  the  burning  rate.  Part  of  this  can 
be  done  experimentally;  for  example,  by  designing  a  bomb  to  burn  a  strand  so  that 
the  method  of  calculating  the  burning  rate  from  the  pressure  time  trace  can  be 
checked.  However,  the  rapid  rate  of  combustion,  as  well  as  the  high  pressures 
and  temperature,  limit  the  amount  of  detailed  knowledge  that  can  be  obtained  from 
experiments.  An  alternative  approach  is  to  attempt  a  mathematical  model  of  the 
ballistic  combustion  cycle*  This  approach  seems  more  feasible  because  the  laws 
of  physics  and  chemistry  are  fairly  well  understood  and  modern  electronic  com¬ 
puters  can  be  used  to  solve  complicated  equations  which  in  the  past  were  unsolv- 
able.  Several  investigations  have  already  been  conducted  and  several  computer 
ballistic  codes  are  now  available.  There  are  the  efforts  of  Kuo1,  Krier  (ref  1), 
Gough  (ref  2),  and  Fisher  and  Trippe  (ref  3).  These  ballistic  codes  emphasize 
the  gas  dynamics  of  the  flame’s  spreading  but  fail  to  consider  the  chemistry 
involved. 


Private  communication,  1976. 
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The  ballistic  codes  are  based  on  an  energy  balance,  the  source  of  which  is 
chemical  in  nature.  The  amount  of  energy  released  at  any  instant  depends  on  the 
energy  of  the  propellant  consumed  minus  the  energy  of  the  products  formed.  The 
products  formed  depend  on  the  pressure,  the  temperature,  and  the  composition  of 
the  propellant  gas.  Thus,  the  temperature,  pressure,  and  propellant  gas  composi¬ 
tion  must  be  accurately  calculated  at  each  step  of  the  simulation.  Calculating 
the  temperature  accurately  requires  knowing  the  heat  capacity  of  the  gas  and  the 
nature  of  the  boundary  layer  which  controls  the  heat  flow  out  of  the  gun.  The 
pressure  must  be  accurately  determined  from  the  temperature  and  gas  composition 
because  it  controls  the  burning  rate  and  performs  the  work  on  the  projectile. 
Thus,  propellant  gas  composition  becomes  a  vital  consideration  in  an  accurate 
model  of  interior  ballistics.  This  composition  is  not  considered  in  the  models 
currently  available . 

Gas  composition  is  also  important  for  other  reasons,  the  most  prominent 
being  muzzle  flash  and  chemical  erosion.  It  is  widely  believed  that  the  most 
significant  contribution  to  gun  erosion  is  thermal;  however,  cool  nitramine  pro¬ 
pellants  have  been  considered  and  found  to  be  more  erosive  than  standard  propel¬ 
lants  with  higher  flame  temperature.  This  could  be  due  to  the  differences  in  the 
transport  properties  of  the  products  or  the  presence  in  the  propellant  gas  mix¬ 
ture  of  reactive  molecules  such  as  H^,  CO,  0^  which  could  be  responsible  for  a 
chemical  attack  on  the  gun  barrel.  Changes  in  formulation  can  be  sought  to  re¬ 
duce  the  erosive  species.  Since  muzzle  flash  is  related  to  the  presence  of  com¬ 
bustible  gases  at  the  muzzle,  knowledge  of  the  gas  composition  and  temperature  at 
each  point  of  the  ballistic  cycle  should  allow  prediction  of  the  extent  of  flash 
for  any  composition. 

Until  now,  gas  composition  has  not  been  a  factor  in  internal  ballistic  cal¬ 
culations  because  it  has  been  difficult  to  include  this  type  of  consideration. 
However,  developments  In  statistical  mechanics  and  in  the  theory  of  simple  fluids 
offer  a  resolution  to  this  problem,  providing  the  products  of  combustion  are 
assumed  to  be  in  thermodynamic  equilibrium.  All  propellants  are  Inherently  un¬ 
stable  substances  and  are  in  metastable  equilibrium  with  their  environment.  The 
combustion  of  a  propellant  Is  not  an  equilibrium  process,  but  If  all  rate  proc¬ 
esses  are  incorporated  into  the  burning  rate  law,  it  may  be  fairly  accurate  to 
assume  that  the  products  of  combustion  are  in  thermodynamic  equilibrium.  This 
can  be  partially  justified  by  estimating  the  number  of  collisions  these  molecules 
undergo  in  a  short  time  compared  to  the  ballistic  cycle.  A  simple  analysis  is 
demonstrated  in  the  appendix.  Turbulence  and  pressure  waves  in  the  gun  tube 
would  also  help  maintain  thermodynamic  equilibrium  between  product  molecules. 

Thus,  to  obtain  realistic  values  for  the  rate  of  propellant  combustion,  it 
is  necessary  to  model  the  combustion  process  precisely.  One  of  the  most  diffi¬ 
cult  aspects  of  this  problem,  even  if  thermodynamic  equilibrium  among  products  is 
assumed,  is  to  obtain  an  accurate  equation  of  state  for  a  propellant  gas  mixture, 
mainly  because  there  is  no  valid  experimental  data  for  a  gas  at  700  MPa  and 
3000 °C.  No  container  has  been  devised  that  will  hold  a  gas  at  this  temperature 
and  at  pressure  long  enough  for  a  static  measurement  to  be  made.  Dynamic  types 
of  measurements  were  made  by  Sandow  (ref  4)  by  burning  propellants  in  a  closed 
bomb.  In  this  experiment  only  the  pressure  was  measured  directly.  A  calculated 
flame  temperature  and  a  constant  moles-of-gas-per-gram  of  propellant  were  used. 
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No  effort  was  made  to  determine  the  extent  of  combustion  or  if  chemical  equilib¬ 
rium  had  been  established.  Consequently,  it  became  necessary  to  calculate  the 
equation  of  state  from  theoretical  considerations.  Fortunately,  the  theory  of 
simple  fluids  also  yields  the  Gibbs  free  energy  of  the  gas  mixture  which  is  nec¬ 
essary  to  calculate  the  equilibrium  constants  between  propellant  gas  molecules 
such  as  the  constant  for  the  water-gas  reaction.  These  constants  are  necessary 
as  a  function  of  pressure  and  temperature  in  order  to  calculate  the  propellant 
gas  composition  at  each  point  in  the  ballistic  cycle. 

This  report  addresses  the  central  problem  of  a  good  equation  of  state  for  a 
propellant  gas.  Some  accurate  experimental  data  (ref  5)  were  obtained  on  N2 
which  is  a  typical  propellant  gas  molecule.  The  data  approach  liquid  densities 
and  a  maximum  temperature  of  1000°C,  the  limit  set  by  conventional  apparatus. 
The  data  are  quoted  to  be  accurate  to  0.3%.  For  an  estimate  of  the  errors  in¬ 
volved  in  using  the  equations  of  state  currently  employed  in  ballistic  calcula¬ 
tions,  this  data  was  fitted  with  empirical  equations  of  state.  Several  theo¬ 
retical  methods  based  on  an  empirical  intermolecular  potential  were  also 
employed.  Then  the  parameters  in  the  empirical  potential  were  adjusted  to  obtain 
the  best  fit  of  experimental  data.  Finally,  the  pressure  of  nitrogen  at  3000  K 
was  calculated. 


COMPUTATIONAL  DETAILS 


Empirical  equations  of  state  were  fitted  to  the  nitrogen  pressure  data  using 
the  non-linear  least— squares  program  of  Moore  and  Zeigler  (ref  6).  This  method 
uses  the  Newton-Raphson  technique  (ref  7)  and  requires  that  derivatives  be  calcu¬ 
lated  with  respect  to  all  adjustable  parameters.  The  Percus-Yevic  equation  (ref 
8)  employs  the  direct  correlation  function  c(r): 

G(r)  =  g(r)-l  =  c(r)  +  p/c( r ' )G( | r-r ' | )dr 
where  g(r)  is  the  radial  distribution  function,  with  the  approximation  that 


c(r)  =  g(r)  {  1  -  exp  [u(r)/k  T] }  . 

u(r)  is  the  intermolecular  potential.  k  and  T  are  defined  in  table  1.  This 
equation  was  solved  using  a  computer  program2  employing  a  method  advocated  first 
by  Broyles  (ref  9).  It  started  with  a  trial  solution  g  and  was  iterated  until  a 
self-consistent  g  was  obtained.  If  the  first  input  value  of  g  is  called  "g^" 
(old)  and  the  result  output  value  is  called  "gout>"  then  the  second  input,  Mg^n" 
(new)  is  given  by: 


The  program  was  supplied  by  Frederic  Mandel,  Univ.  of  Cincinnati  College  of 
Medicine,  Dept,  of  Pharmacology  and  Cell  Biophysics,  Cincinnati,  Ohio. 
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gln(new)  “  a  gln(old)  +  (1-a  )  gQut 

with  0<a  <1.  At  high  densities  it  is  often  necessary  to  use  large  values  of  a  to 
achieve  convergence . 

The  third  virial  coefficients  were  extrapolated  from  the  tables  compiled  by 
Bird  et  al  (ref.  10).  Attempts  to  compute  these  values  using  modern  methods  of 
evaluating  multidimensional  integrals  of  Haselgrove  (ref  11)  and  Conroy  (ref  12), 
were  unsuccessful  because  of  inaccuracies  in  the  methods.  A  computer  program  to 
perform  the  Verlet-Weis  perturbation  theory  (ref  13)  was  then  written.  A  least- 
squares  computer  program  which  does  not  use  derivatives  was  used  in  order  to 
adjust  the  Lennard-Jones  potential  for  the  best  fit  of  nitrogen  pressure  data. 
This  technique  was  programmed  incorporating  Powells  ideas  (ref  14)  concerning 
the  variation  of  parameters  along  conjugate  directions. 


RESULTS  AND  DISCUSSION 

In  the  investigation  of  the  role  of  the  equation  of  state  in  propellant 
combustion,  the  pressure-volume-temperature  data  of  Malbrunot  and  Vodar  (ref  5) 
was  used.  These  data,  stated  to  be  accurate  to  0.3%,  span  the  temperature  range 
of  200°C  to  1000°C  and  0.1618  g/cra3  to  0.8222  g/cm3  in  density.  These  data  were 
fit  to  popularly-used  equations  of  state  at  various  temperatures  to  ascertain  the 
degree  of  accuracy  of  the  equations  of  state  currently  employed  in  interior  bal¬ 
listic  calculations.  The  data  at  600°C  were  excluded  because  of  obvious 
errors.  Consideration  was  given  to  the  virial  expansion  using  just  the  second 
virial  coefficient  as  employed  in  the  Calspan  code  (ref  3),  the  Able-Noble  equa¬ 
tion  which  is  used  by  other  ballisticians,  and  two  virial  coefficients  as  used  in 
the  BLAKE  code  (ref  15).  Also  tried  was  an  equation  of  state  which  had  recently 
been  found  to  be  very  accurate  for  hard  spheres  (ref  16).  The  functional  form  of 
the  equations  is  summarized  in  table  1.  Tables  2  through  9  record  the  calculated 
results,  the  experimental  values  of  Malbrunot  and  Vodar,  and  the  errors  reported 
as  a  percentage  of  the  experimental  values.  It  was  evident  that  the  need  exists 
for  a  simple,  accurate  equation  of  state  that  can  be  used  in  ballistic  calcula¬ 
tions  as  well  as  a  method  which  will  becomes  more  accurate  as  temperatures  in¬ 
crease.  Since  there  is  no  experimental  data  at  3000  K,  a  search  should  be  made 
for  a  simple  method  for  extrapolating  to  high  temperature. 

Calculations  in  tables  2  through  9  showed  that  none  of  the  simple  methods 
works  well.  The  Able-Noble  equation  was  particularly  disappointing.  This  ap¬ 
proach  worked  at  200°C,  but  as  temperature  increased,  agreement  with  the  exper¬ 
iment  steadily  became  worse,  until  at  1000°C  there  were  9%  errors.  This  is  under 
the  best  situation  where  the  covolume  was  allowed  to  be  a  function  of  temperature 
(covolume  increased  with  temperature).  Usually  the  covolume  is  assumed  to  be  a 
constant,  and  in  this  case  errors  as  great  as  15%  can  be  encountered.  In  a  real¬ 
istic  computer  model,  errors  of  this  magnitude  cannot  be  tolerated.  To  obtain 
reasonable  results,  it  is  necessary  to  adjust  "shot-start"  pressure,  burning  rate 
parameters,  and  heat  losses. 

Keeping  only  the  "second"  virial  coefficient  is  also  a  very  poor  approxima¬ 
tion  although  this  approach  improves  with  temperature.  Using  the  "second"  and 
"third"  virial  coefficients  gives  the  best  results,  but  large  increases  of  the 
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"third"  virial  coefficient  with  temperature  make  any  attempt  at  extrapolation  to 
3000  K  questionable.  The  use  of  the  equation  of  state  which  is  known  to  be  accu¬ 
rate  for  hard  spheres  seems  to  yield  readily  to  physical  interpretation.  This 
equation  gives  poor  results  at  200°G  where  the  attractive  part  of  the  potential 
is  still  important.  As  temperature  increases,  these  attractive  forces  are  more 
easily  overwhelmed  and  the  molecules  start  to  behave  like  hard  spheres.  Also, 
the  effective  hard-sphere  volume  decreases  as  a  function  of  temperature.  When 
plotted  as  a  function  of  the  reciprocal  of  the  absolute  temperature,  it  extrap¬ 
olates  smoothly  to  a  value  which  can  be  used  at  3000  K.  This  extrapolation  is 
depicted  in  figure  1. 

Because  of  the  rather  poor  performance  of  the  empirical  equations  of  state, 
the  theories  of  simple  fluids  were  applied  to  calculate  the  pressure  of  nitrogen 
gas  and  compare  the  calculated  values  to  the  experimental  data.  These  methods 
employ  an  intermolecular  potential  function.  The  Percus-Yevic  method  calculates 
thermodynamic  properties  in  terms  of  the  two-particle  distribution  function.  The 
perturbation  theory  of  Verlet-Weis  and  the  more  conventional  virial  expansion 
were  used  to  calculate  thermodynamic  properties  directly.  The  virial  expansion  is 
exact  and  becomes  approximate  because  only  the  second  and  third  virial  coeffi¬ 
cients  were  used.  Results  of  computations  employing  the  Lennard-Jones  potential 


are  shown  in  tables  10  through  17.  In  this  study, 

a  =  0.3698  nm  and  r  =  95.05  K  for  N„ . 

k  2 

Several  conclusions  can  be  made  on  the  basis  of  these  calculations.  The  most 
accurate  pressures  were  obtained  with  the  Percus-Yevic  equation.  The  greatest 
deviation  was  7.6%  at  low  temperature  and  high  density.  The  Verlet-Weis  theory 
was  less  accurate.  Since  all  the  methods  improve  with  increasing  temperature, 
the  pressure  of  nitrogen  was  computed  at  3000  K.  The  results  of  these  calcula¬ 
tions  are  recorded  in  table  18  and  compared  with  extrapolated  values  obtained 
using  the  hard-sphere  equation  of  state.  Good  agreement  was  obtained  with  the 
Percus-Yevic  equation.  The  agreement  was  reduced  when  the  Verlet-Weis  perturba¬ 
tion  theory  was  employed.  (The  two-term  virial  equation  is  always  too  low  be¬ 
cause  higher  order  virial  coefficients  are  not  used.) 

For  accurate  interior  ballistic  calculations,  agreement  of  several  percent 
is  inadequate.  Good  agreement  between  the  Percus-Yevic  and  the  Verlet-Weis  the¬ 
ories  indicated  that  the  discrepancies  were  probably  due  to  the  fact  that  the 
Lennard-Jones  potential  did  not  accurately  describe  the  nitrogen  intermolecular 
potential.  Since  the  o  and  the  e  are  numbers  obtained  by  fitting  experimental 
data,  these  quantities  were  adjusted  to  obtain  the  best  fit  of  the  experimental 
nitrogen  pressures.  The  Percus-Yevic  theory  is  computationally  too  complex  for 
this  type  of  effort,  but  the  Verlet-Weis  theory  requires  little  computer  time  and 
is  easily  adapted  to  this  type  of  approach.  The  deviation  between  experimental 
and  calculated  points  is  some  function  of  o  and  e;  however,  this  function  is 
complex  and  there  does  not  seem  to  be  any  simple  way  to  calculate  derivatives 
with  respect  to  a  and  e .  For  this  reason,  the  Powell  method  was  used  which  re¬ 
quires  no  derivatives  to  minimize  the  deviation  with  respect  to  a  and  e.  The 
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best  values  of  a  and  e  were  obtained  for  each  temperature;  results  are  shown  in 
table  19. 


The  calculations  fall  into  two  categories.  In  the  temperature  range  from 
2Q0°C  to  600°C,  good  agreement  is  obtained  and  the  distortion  of  the  potential  is 
small.  At  elevated  temperatures,  the  agreement  is  poor  and  distortion  of  the 
Lennard-Jones  potential  is  much  greater.  This  illustrates  a  well-known  inade¬ 
quacy  of  the  Lennard-Jones  potential;  i.e.,  the  repulsive  part  is  too  steep.  At 
high  temperatures,  the  attractive  part  is  not  very  important  and  the  e  can  be 
adjusted  to  obtain  a  more  realistic  repulsive  potential,  although  not  very  suc¬ 
cessfully.  To  remedy  this  problem  a  more  realistic  potential  is  required.  A 
modified  version  of  the  Lennard-Jones  6-12  potential  has  been  suggested  by 
Maitland  and  Smith  (ref  17): 

“<**>  -  e  1  -  (jp*)  Cr*r6l, 

m  =  13  +  Y  (r*-l) , 

where  r*  =  r/rm^n  and  rm^n  is  the  separation  at  the  minimum  energy,  -e .  This 
potential  is  of  the  Lennard-Jones  6-12  type  if  m  is  set  equal  to  12.  A  more 
flexible  and  realistic  potential  is  achieved  by  allowing  the  m  to  vary  with  in- 
termolecular  distance.  Aziz  (ref  18)  found  that  if  this  potential  is  applied  to 
neon,  the  parameters  e  ,  rmin,  and  y  can  be  evaluated  from  three  dilute  gas  bulk 
properties;  thermal  diffusion,  viscosity,  and  second  virial  coefficients.  The 
resulting  potential  predicts  not  only  dilute  gas  bulk  properties  but  also  spec¬ 
troscopic  and  molecular  beam  data  as  well  as  much  more  complex  realistic  poten¬ 
tials  (ref  17).  This  new  potential  could  yield  accurate  pressure  determinations 
for  propellant  gases  at  high  temperature. 


CONCLUSIONS 


The  simple  equations  of  state  usually  employed  in  ballistic-type  calcula¬ 
tions  such  as  the  Able-Noble  equation  and  the  two-term  virial  expansion  are  not 
accurate.  The  virial  expansion  is  better  for  ballistic  calculations  since  it 
improves  with  temperature,  but  it  is  not  amenable  to  accurate  extrapolations.  On 
the  other  hand,  modern  theories  of  fluids  such  as  the  Percus-Yevic  theory  and  the 
perturbation  theory  of  Verlet-Weis  offer  a  reasonable  method  of  determining  an 
equation  of  state  for  propellant  gases  at  high  temperature  and  pressure,  with  the 
best  results  being  consistently  obtained  with  the  Percus-Yevic  equation.  The 
Lennard-Jones  intermolecular  potential  is  not  flexible  enough  for  description  of 
nitrogen-nitrogen  interactions  at  high  temperatures. 
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RECOMMENDATIONS 


Because  of  the  poor  performance  of  the  equations  of  state  used  in  convention¬ 
al  ballistic  codes  and  the  total  lack  of  consideration  of  the  propellant  gas 
composition,  there  may  be  merit  in  combining  the  BLAKE  thermodynamic  code  with  an 
interior  ballistic  code.  This  combination  would  permit  calculation  of  the  gas 
composition  at  each  point  in  the  ballistic  cycle  and  would  provide  an  equation  of 
state  based  on  a  two-term  virial  expansion  which  is  more  accurate  than  the  equa¬ 
tion  of  state  presently  employed.  Such  a  program  could  be  combined  with  a  least- 
squares  technique  to  obtain  burning  rate  data  from  a  closed  bomb  pressure  time 
trace,  the  burning  rate  parameters  being  adjusted  until  the  best  least-squares 
fit  of  the  pressure  time  trace  is  obtained. 

The  new  molecular  interaction  potential  of  Maitland  and  Smith  should  be 
utilized  to  get  a  more  accurate  pressure  calculation  for  nitrogen  as  well  as 
other  propellant  gases.  This  could  result  in  calculations  of  more  accurate 
thermodynamic  properties  of  propellant  gas  mixtures. 

The  experimental  data  obtained  by  Malbrunot  and  Vodar  for  nitrogen  did  not 
include  pressure  information  above  1000°C,  but  it  proved  helpful  in  this  study. 
Similar  experimental  pressure  data  should  be  obtained,  collectively  and  sepa¬ 
rately,  on  other  propellant  gas  molecules  such  as  CO,  H^,  H20,  C02,  N2,  and  CH^. 
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APPENDIX 

PARTIAL  JUSTIFICATION  OF  CHEMICAL  EQUILIBRIUM  OF  PROPELLANT  GASES 
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For  achievement  of  chemical  equilibrium  in  a  propellant  gas  during  the  bal¬ 
listic  event,  several  criteria  must  be  met.  There  should  be  sufficient  colli¬ 
sions  to  give  individual  molecules  a  chance  to  react  with  each  other  in  a  short 
time  compared  to  the  total  ballistic  cycle.  The  temperature  should  be  high 
enough  for  chemical  reaction;  that  is,  the  temperature  should  be  above  the 
"freeze-out"  temperature.  Since  at  high  density,  mean-free  paths  of  molecules 
may  be  quite  short,  mixing  of  propellant  gas  molecules  may  be  very  slow.  It  is 
therefore  vital  that  the  contents  of  the  chamber  be  inherently  homogeneous. 

For  estimating  the  number  of  collisions  in  a  mixture  of  propellant  gases, 
formula  for  the  number  of  collisions,  f  was  used: 

2  . 

f  '  i  d  n  /2  v 

where  d  is  the  diameter  of  the  molecule,  v  is  its  velocity  and  n  is  the  number 
density.  For  a  typical  molecule  at  room  temperature,  d  is  2.25  x  10-8  cm,  and 

the  velocity  is  10s  cm/sec.  At  100  KPa  of  pressure,  there  are  about  3  x  1019 

molecules/cm3.  Substitution  of  these  values  in  the  above  formula  resulted  in  7  x 
109  collisions  per  second.  If  the  time  for  the  ballistic  cycle  is  taken  as  10 
milliseconds  and  the  total  is  divided  into  109  steps  for  the  purpose  of  numer¬ 
ical  integration,  in  each  of  the  small-time  intervals  there  are  7  x  103  colli¬ 

sions.  There  should  be  sufficient  collisions  to  maintain  chemical  equilibrium. 
At  increased  temperature  and  pressure  (as  is  actually  experienced  in  the  gun 
environment),  both  v  and  n  increase  so  that  even  more  collisions  occur  to  main¬ 
tain  chemical  equilibrium. 

Freeze-out  temperature  for  a  propellant  is  about  1000  K  while  the  flame 
temperature  is  over  2000  K,  so  that  there  should  be  sufficient  energy  in  an  aver¬ 
age  collision  to  insure  the  possibility  of  chemical  reaction  and  the  maintenance 
of  chemical  equilibrium.  Because  mean-free  paths  of  propellant  gas  molecules  may 
be  quite  short  at  high  densities,  mixing  may  be  very  slow  between  different  re¬ 
gions.  However,  if  the  propellant  is  homogeneous  and  the  products  of  combustion 
are  the  same  in  all  parts  of  the  chamber,  extensive  mixing  should  not  be  neces¬ 
sary.  Any  turbulence  of  pressure  waves  that  travel  in  the  combustion  chamber 
will,  of  course,  promote  mixing  if  it  is  required. 
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Table  1.  Empirical  and  theoretical  equations  of  state 
currently  employed  in  propellant  calculations 


Able-Noble 

p  =  pk  T(l-y)- 1 


Virial  Expansion* 

p  =  pkT(l  +  B2(T)  p  +  B3(T)p2  +  - 

Hard  Sphere 


p  =■  pkT(l  +  y  +  y2  -  y3)(l  -  y)  3 
p  =  gas  density  (g/cm3) 

b  =  adjustable  parameter  (cm3/g) 

y  -  Pb 

k  ■  Baltzman  constant  (MPa  cm3/g  K) 

P  =  pressure  (MPa) 

T  =  temperature  (K) 


*  This  equation  can  be  either  empirical  or  theoretical  depending  on  whether 
B(T),s  are  determined  from  experimental  data  or  calculated  from  an  interraolecu- 
lar  potential  as  employed  in  BLAKE. 


14 


Table  2.  Observed  and  calculated  pressure  of  nitrogen 
at  200°C  using  empirical  equations  of  state 
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Table  3.  Observed  and  calculated  pressure  of  nitrogen 
at  300°C  using  empirical  equations  of  state 
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Table  5.  Observed  and  calculated  pressure  of  nitrogen 
at  500°C  using  empirical  equations  of  state 
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Table  6.  Observed  and  calculated  pressure  of  nitrogen 
at  700°C  using  empirical  equations  of  state 
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Table  7.  Observed  and  calculated  pressure  of  nitrogen 
at  800°C  using  empirical  equations  of  state 
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Table  8.  Observed  and  calculated  pressure  of  nitrogen 
at  900°C  using  empirical  equations  of  state 
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Table  9.  Observed  and  calculated  pressure  of  nitrogen 
at  1000°C  using  empirical  equations  of  state 
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Table  10.  Observed  and  calculated  pressure  of  nitrogen  at  200 °C 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Table  11.  Observed  and  calculated  pressure  of  nitrogen  at  300°C 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Properties  of  Nitrogen  up  to  1000°C  and  5000  Bar,”  Physica,  vol  66,  1973,  p  351. 


Table  12.  Observed  and  calculated  pressure  of  nitrogen  at  400°C 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Experimental  data  taken  from  P.  Malbrunot  and  B.  Vodar,  "Experimental  pVT  Data  and  Thermodynamic 
Properties  of  Nitrogen  up  to  1000°C  and  5000  Bar,"  Physica,  vol  66,  1973,  p  351. 


Table  13.  Observed  and  calculated  pressure  of  nitrogen  at  500°C 
using  the  Lennard-Jones  6-12  intermole cular  potential 
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Table  14.  Observed  and  calculated  pressure  of  nitrogen  at  700 °C 
using  the  Lennard-Jones  6-12  interraolecular  potential 
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Table  15.  Observed  and  calculated  pressure  of  nitrogen  at  800°C 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Table  16.  Observed  and  calculated  pressure  of  nitrogen  at  900°C 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Table  17.  Observed  and  calculated  pressure  of  nitrogen  at  1000°' 
using  the  Lennard-Jones  6-12  intermolecular  potential 
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Table  18.  Calculation  of  nitrogen  pressure  at  T  =  3000  K 
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Table  19.  Fit  of  Lennard-Jones  6-12  potential  to  nitrogen 
pressure  data  using  Verlet-Weis  theory 
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